Phase-Locking Transition in a Chirped Superconducting Josephson Resonator 



00 

o 
o 

(N 

>>: 



CD 



o 

(N 
> 

ON 

in 

ON 

in 
o 

00 

O 



X 



O. Naaman,^ J. Aumentado,^ L. Friedland,'^ J. S. Wurtele,^'^ and I. Siddiqi^ 

^ Quantum Nanoelectronics Laboratory, Department of Physics, 
University of California, Berkeley, California 94720, USA 
^National Institute of Standards and Technology, 325 Broadway, Boulder, Colorado 80305, USA 
'^Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel 
"^Department of Physics, University of California, Berkeley, California 94720, USA 
^Center for Beam Physics, Lawrence Berkeley National Laboratory, Berkeley CA 94720, USA 

(Dated: May 14, 2008) 



PACS numbers: 85.25.Cp, 05.45.-a, 03.67.Lx, 74.50.+r 



By coupling a harmonic oscillator to a quan- 
tum system it is possible to perform a dispersive 
measurement^'^ that is quantum non-demolition 
(QND), with minimal backaction^. A non-linear 
oscillator has the advantage of measurement gain, 
but what is the backaction? Experiments on su- 
perconducting quantum bits (qubits) coupled to 
a non-linear Josephson oscillatoii^>^i^ have thus 
far utilized the switching of the oscillator near a 
dynamical bifurcation for sensitivity'"'^, and have 
demonstrated partial QND measurement. The 
detailed backaction associated with the switching 
process is complex, and may ultimately limit the 
degree to which such a measurement can be QND. 
Here we demonstrate a new dynamical effect in 
Josephson oscillators by which the bifurcation can 
be accessed without switching. When energized 
with a frequency chirped drive with an amplitude 
close to a sharp phase-locking threshold, the os- 
cillator evolves smoothly in one of two diverging 
trajectories^ — a pointer for the state of a qubit. 
The observed critical behavior agrees well with 
theory and suggests a new modality for quantum 
state measurement. 

A superconductor-insulator-superconductor tunnel 
junction is a unique electrical circuit element that can 
be simultaneously non-linear and non-dissipative. The 
Josephson equations parameterize the non-linear tunnel 
current / = /osinJ and voltage Vj — hS/2e in terms 
of the gauge-invariant superconducting phase difference 
6 across the junction and its time derivative S. These 
equations describe a non-linear inductor with inductance 
Lj{S) — h/{2eIoCOsS), which can be shunted with 
a low-loss reactanceiS to form a high quality-factor 
(Q) anharmonic oscillator. Coupling this Josephson 
oscillator to a qubit produces a state dependent shift 
of the resonant frequency — thus realizing a dispersive 
measuremen t 1 1 1 1 ^ 1 13 ^ When probed with a small num- 
ber of photons, the Josephson oscillator is essentially 
harmonic and the measurement does not project the 
quantum state out of the qubit basis: it is therefore 
quantum non-demolition (QND), and has a well charac- 
terized minimal backactiort^'^ i^'^'^^ . When driven more 
strongly, anharmonicity causes the resonant frequency 



to vary with the amplitude of oscillatior*i^, enhancing 
the oscillator's driven response to a frequency shift and 
thus "amplifying" the qubit "signal" . Since the junction 
has vanishing internal dissipation, this process may also 
result in minimal decoherence. Moreover, the Josephson 
oscillator can also bifurcate^ into two metastable os- 
cillation states — one which is low-amplitude and phase 
lagging and one which is high-amplitude and phase 
leading — resulting in a projective measurement. This 
has been demonstrated in the Josephson bifurcation 
amplifier (JBA)^ in which the oscillator is probed with 
a fixed frequency drive with varying amplitude^i^. The 
dynamical state of the oscillator depends on the qubit 
state, and measurement sensitivity arises from a sharp 
threshold for switching between the two dynamical 
states. While such a non-linear oscillator has the 
advantage of measurement gain, the precise nature of 
its backaction on the quantum system, especially that 
associated with the switching process, and the degree 
to which a QND measurement is possible^ is still being 
investigated. 

We demonstrate a new measurement scheme that still 
drives the Josephson oscillator to the bifurcation regime 
but does not involve any switching process. Instead 
of applying an amplitude modulated drive, we apply 
a chirped, microwave frequency drive to access a phe- 
nomenon known as autoresonancei^. In response to a fre- 
quency modulation, the oscillator may either phase-lock 
to the chirped carrier and latch to the high-amplitude 
oscillation state, or not lock to the drive and remain in 
the small-oscillation state. These two outcomes are sep- 
arated in parameter space by a sharp thresholdi^i that 
scales with the chirp rate and is sensitive to the junction 
Iq. Throughout its evolution, the system tracks a single 
basin of attraction — no switching occurs, thus avoiding 
any potential backaction associated with switching dy- 
namics and transient oscillations. We call this new mea- 
surement device a Josephson chirped amplifier (JCA). 
We show that the observed threshold behavior is in excel- 
lent agreement with both analytic theory and numerical 
simulations, and estimate the performance of the JCA as 
a quantum state readout device. 

In our experiment, the Josephson oscillator was formed 
by an Al/AlO^j/Al Josephson tunnel junction (Fig. [TJ;) 
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FIG. 1: Experimental setup and circuit model, a, 

Schematic of the experimental setup, b, Series RLC model of 
the resonator, c, Scanning electron micrograph of the junc- 
tion. 



placed in the middle of a high-Q Nb half- wave copla- 
nar waveguide resonator with a characteristic impedance 
Zo = 50 ri, and symmetrically coupled to the 50 en- 
vironment via capacitors Cc ~ 10 fF. The measurements 
were performed in a dilution refrigerator at T=20 mK; 
the experimental setup is shown schematically in Fig. 
[T^. Microwave excitation was applied to the resonator 
using an HP8780A vector signal generator, frequency- 
modulated with a triangle waveform to provide a linear, 
phase-continuous frequency chirp with a 50 MHz span. 
Typical chirp rates, a = duj/dt, ranged from a/27r = 10^^ 
Hz/s to 10^3 Hz/s. The transmitted microwave signal, 
Vout, was amplified and then demodulated (1.9 MHz IF 
bandwidth) to find its amplitude and phase. The ampli- 
tude of current oscillations in the resonator is given by 
(2|Kut|/Zo) X v/QA. 
We first measure the oscillator in steady state (Fig. 
[2^). From the microwave transmission. Pout) at low 
power in the linear regime (Fig. [2^, blue curve), we mea- 
sure ijJo/27r = 1.61564 GHz and infer a quality factor 
Q = 27500 ± 1000. The Q is limited by coupling to 
the 50 fl environment and is within 10% of its predicted 
value Q = tt / AZqUJqC'^ . When the drive power exceeds 
a critical value Pc — —148 dBm, the oscillator response 
bifurcates into two branches^. From the measured Pc we 
estimate the junction critical current Iq — 0.61±0.04 ^A; 
this value is consistent with room temperature dc re- 
sistance measurements on co-fabricated junctions. The 
stated uncertainties are dominated by possible cryogenic 
variations in the attenuation of the coaxial lines attached 
to the resonator. Additionally, the observed power at 
which bifurcation occurs at different drive frequencies 
agrees well with theoryiS using the measured values of 
Pc and Q. 

Having characterized its parameters, we proceed to ex- 
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FIG. 2: Microwave transmission and phase space tra- 
jectories, a, Transmitted power Pout in steady state, mea- 
sured with a network analyzer. Input powers are —155.5 dBm 
(blue) to —123 dBm at 2.5 dB step. The system exhibits bi- 
furcation above a critical power Pc = —148 dBm (1.58 x 10~^* 
W). b, Simulated in-phase and quadrature components of cur- 
rent oscillations for detuning {Lj—ujQ)/2n = 0.5 MHz and drive 
power ramped in time from -126 dBm to -120 dBm. A and B 
label the low- and high- amplitude attractors, 8 is the phase 
of the oscillations, c, Response of the resonator to a chirped 
drive, a/2TT = 5 x lO" Hz/s, at Pi„ = -126 dBm (blue) and 
— 120 dBm (red), d, Simulated quadratures of the current in 
a chirped resonator with the same parameters as the data in 
c. The simulation was terminated at ui/2tt = 1.6106 GHz. 



cite the oscillator with a chirped microwave drive. Fig.[2j; 
shows the resonator response to a downwards frequency 
chirp at the rate of a/27r = 5 x 10^^ Hz/s. At low 
power (blue) the amplitude of oscillations initially in- 
creases as the chirp passes through ujq. However, as the 
chirp progresses towards lower frequencies, the resonator 
decouples from the drive and rings down to rest. For 
a stronger drive amplitude (red) the response of the res- 
onator changes dramatically: as the chirp passes through 
ujo, the resonator phase becomes locked to the drive and 
its amplitude grows with time. A threshold for phase- 
locking can be seen in Fig. [3l which shows the normal- 
ized amplitude of current oscillations, I/Iq, as a function 
of frequency and drive power for a fixed chirp rate. In 
region 1, at low powers, no phase locking is observed 
and the junction current grows only in the vicinity of 
luq- Above a critical drive power, Pc in Fig. [Sj the res- 
onator remains phase locked and its amplitude continues 
to grow up to a deterministic maximum, set either by 
damping (region 2) or by Iq (region 3). The existence 
of a threshold for phase-locking, which obeys a univer- 
sal scaling law, was first observed in the context of non- 
neutral plasmas^; our work is the first observation of this 
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FIG. 3: Microwave response of a chirped non-linear 
resonator. Normalized amplitude of current oscillations in 
the resonator as a function of drive power and frequency, with 
the frequency swept downwards at a rate of 5 x 10^^ Hz/s. 
The arrow indicates the critical power for phase-locking. 
In region 1 where Pin is sub-critical, the oscillation amplitude 
grows in the vicinity of uo and rings back to zero below ujq. 
Above Pc, the oscillator phase- locks and its amplitude grows 
until it reaches a maximum set either by damping (region 2) 
or by Jo (region 3). 

transition in a micro-electronic circuit operating at GHz 
frequency and mK temperature — five orders of magni- 
tude higher in frequency and lower in temperature than 
Ref. i. 

We proceed to briefly analyze the threshold phe- 
nomenon (a detailed analysis will be given elsewhere). 
The dynamics of a Josephson junction, ac biased through 
a resonant cavity near its resonance, can be modeled by 
the equivalent series RLC circuit shown in Fig. [TJaiSii^, 
where L = L + Lj(0) w ttZq/2loq = 7.74 nH, C « 
II-kZqujq = 1.25 pF, and cjg = l/LC. The effective 
loss due to external loading is i?off = ttZq/2Q, and 
the effective drive amplitude is Vd = Vin\/'n'/Q where 
Vin — ^/2ZoPin is the amplitude of the incident wave 
referred to a matched load. From Kirchoff's voltage 
law we obtain an equation of motion for the charge 
q{t) — J I{t')dt' in the circuit. Defining the non-linearity 
ratio 2a — Lj{Q)/L, and the dimensionless time r = ujot, 
charge q — {ujQ\/2a/ Ii^)q^ and drive e = (•\/2a/Z/oa-'o)Vd, 
and expanding the voltage drop across the junction to 
second order in ///q, the equation of motion in the 
weakly non-linear regime^, becomes: 

g(l + gV2) + = ecos0d, (1) 

where v = l/Q, (fid = — at^ /2 is the phase of the 
chirped drive, and the derivatives are with respect to t. 

Transforming to a chirped frame rotating with the 
drive [where equation ([1]) reduces to the Duffing model^"], 
and neglecting fast oscillating terms, we cast equation ([T]) 



into an equation for the complex variable 4* = Aexp(i0), 
where A is proportional to the amplitude of ^(t), and 
is the phase mismatch between q{T) and the drive: 

+ + i7/2)* = /i- (2) 

Here f = a-^/^r/cjo, 7 ~ ^qoT^^"^ IQ, and = 
t u?^"^ a^'^l^ 1%. Equation ([2]), with the initial condition 
^ = Oatf-^-— oo, admits two asymptotic solutions 
at f — > -|-oo: a solution decoupled from the drive for 
/i < /icr, in which grows quadratically with time, and 
a phase locked solution = for \i > [icr- 

In the absence of damping^" the transition between the 
two solutions of equation ([2]) occurs for = 0.41. When 
damping is present, increases and becomes dependent 
on the damping: ^icril) — Mcr(l + 07 -|- 67^) for 7 ^ 1, 
with coefficients a — 1.06 and b = 0.67 found numerically. 

Combining the definitions above for fi and e, we find 
the critical drive amplitude for phase locking: 




The scaling law V^"^ oc a^^^ is exact only in a lossless res- 
onator. For finite dissipation, 7 cx a~^^'^ and corrections 
to this law enter through fieri"/) ■ 

To compare the observed value of the phase-locking 
threshold, T/f, to the theory outlined above, we first 
note that the probability of phase-locking, Piock, ob- 
tained from averaging the response from 5,000 frequency 
sweeps^^, grows from zero to one over a finite range of 
drive amplitudes. A typical phase-locking probability 's- 
curve' is plotted in the inset of Fig. 2] as a function of Vd, 
for a/27r = 2 x 10^^ Hz/s, where we define the threshold 
at Piock = 0.5. We measured V^^ as a function of chirp 
rate, and the results are plotted in Fig. |4] (□). We also 
performed numerical simulations of the fully non-linear 
equation of motion for q{t) in the series RLC model of 
Fig. [T]3, without making any approximations to the volt- 
age drop across the junction, and using experimentally 
determined parameters. We simulated both a resonator 
with Q = 2.75 x 10^ and a lossless one. The numerical 
simulations (-I-) and the predictions from equation ([3]) 
with (solid line) and without (dashed line) damping are 
also plotted in Fig. [H We observe excellent agreement 
between experiment and theory if damping is included. 
Note that this is not a fit — all parameters are fixed to 
their experimentally determined values. Moreover, the 
agreement of equation ([3]) with the fully non-linear sim- 
ulations confirms the validity of the weak non-linearity 
assumptions made in the threshold analysis — the phase- 
locking transition occurs at small amplitudes. 

From the measured width of the phase-locking thresh- 
old we can estimate the potential sensitivity of the JCA, 
where as a benchmark we consider detecting a 1 % vari- 
ation in Jq, a typical signal associated with the transi- 
tion between the ground and excited state of a super- 
conducting 'Quantronium' qubit^'^. The discrimination 
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FIG. 4: Critical drive scaling and phase-locking prob- 
ability curve. Comparison of the critical drive voltage Vd 
obtained from the experiment (□) to equation ((3| with (solid 
line) and without (dashed line) dissipation. Also shown are 
values for the critical drive from simulations of the fully non- 
linear equation of motion for the equivalent RLC circuit (+), 
for lossless resonator (blue) and for Q = 2.75 x 10* (red). 
Inset: probability of phase-locking for drive amplitudes near 
threshold at a/2'K = 2 x 10^^ Hz/s, evaluated over a 0.5 MHz 
frequency band centered at 1.6095 GHz. Solid line is a sig- 
moidal fit. 

power of the device can be estimated from the fractional 
change in critical current, AIq/Iq, that will shift the 
's-curve' by an amount equal to the threshold width: 
A/o/Jo = 2AV^''/3V^'' from equation ©. For the data 
in the inset of Fig. [4l and defining the width AT^'^'' for 
0.27 < Piock < 0.73, we find A/q/Zo = 9.6 x IQ-^. The 
1 % variation in Iq (6.1 nA for our device) can be re- 
solved with ~ 46 % contrast; the ultimate sensitivity of 
the JCA requires an understanding of the dependence of 
the threshold width on noise, both classical and quan- 
tum, which is currently being pursued. For this chirp 
rate of a/2Tr = 2 x 10^^ Hz/s, a single measurement can 
be accomplished in less than 10 fis. Significantly faster 
chirp rates and shorter measurement times are in princi- 
ple possible. 

The dynamics associated with operating the amplifier 
with a frequency modulated (FM) drive versus an am- 
plitude modulated (AM) drive are quite different. In 
Fig. ^p,d we show the calculated in-phase and quadra- 
ture phase components of I/Iq for AM and FM drives, 
respectively. In the AM case, for detuning uj < ujq, the 
system can switch from a low-amplitude oscillation at- 



tractor (A) to a high-amplitude oscillation one (B) , with 
significant phase oscillations associated with the switch- 
ing as shown in Fig. [^b. In fact, for large enough de- 
tuning, the junction's phase can exceed 6 = tt, which 
can induce state mixing when measuring a quantum sys- 
tem. Trajectories for the chirped drive are shown in Fig. 
[5Jl, where the two traces correspond to the experimental 
data shown in Fig. [21;. For sub-critical drives (blue), the 
system builds amplitude as the chirp frequency crosses 
loq, but remains in the low-oscillation state and eventu- 
ally relaxing back to zero amplitude. For super-critical 
drive (red), the system smoothly follows the evolution 
of the high-amplitude attractor without spurious oscilla- 
tions of the phase. Another feature of the JCA is that 
we can potentially perform quantum measurement with 
very few photons (near loq), but "latch" and record the 
signal with a large number of photons at a final frequency 
Lu < luq, so that minimal fidelity is lost due to noise in 
the measurement electronics. 

In summary, we have observed that a high-Q res- 
onator embedding a Josephson junction subject to a 
chirped drive exhibits a transition to phase-locking with 
a sharp threshold that can be used for detecting small 
changes in Iq. This device — the Josephson chirped 
amplifier — provides a new route for non-linear, disper- 
sive quantum measurement and a test bed for non- 
equilibrium quantum statistical mechanics in a chirped 
reference frame. We studied the dependence of the criti- 
cal drive on the chirp rate and found excellent agreement 
with both theory and simulations. The phase- locked 
highly excited state and the associated phase-locking 
threshold are a robust feature of chirped non-linear os- 
cillators, and should be observable in any system with 
low loss and weak non-linearity: electrical, mechanical, 
or photonic, for example. Finally, we note that another 
approach to observing these phase-locking effects is to 
modulate the resonant frequency of the oscillator in the 
presence of a harmonic drive; a candidate technology for 
such tuning was recently demonstrated in Ref. t23j in the 
context of variable coupling elements for qubits. 
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